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Abstract 

Vacancy-induced magnetism in graphene bilayers is investigated using spin-polarized density 
functional theory calculations. One of two graphene layers has a monovacancy. Two atomic 
configurations for bilayers are considered with respect to the position of the monovacancy. We 
find that spin magnetic moments localized at the vacancy site decrease by 10% for our two 
configurations, compared with the graphene monolayer with a monovacancy. The reduction of the 
spin magnetic moment in the graphene bilayers is attributed to the interlayer charge transfer from 
the adjacent layer to the layer with the monovacancy, compensating for spin magnetic moments 
originating from quasilocalized defect states. 
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I. INTRODUCTION 



Vacancy defects in crystalline solids have been of fundamental interest in materials science 
and condensed matter physics. Of particular interest, both technologically and academically, 
are the vacancy defects in nanostructures of sp^-bonded carbon based nanostructures [1,2]. 
They exhibit room-temperature (RT) ferromagnetism [3], and the underlying physics of their 
ferromagnctism is different from that of conventional ferromagnetic metals such as iron and 
cobalt [4,5]. Successive progress in the understanding of ferromagnetism in carbon-based 
nanostructures has been achieved by theoretical approaches based on spin-polarized density 
functional theory [6-9] and experiments on irradiated carbon nanostructures [10-12]. To il- 
lustrate, for a graphene monolayer (a two-dimensional honeycomb lattice of C atoms) with a 
monovacancy, it has been elucidated that spin magnetic moment is localized at the vacancy 
site and that RT-fcrromagnetism in the graphene monolayer originates from localized sp^ 
dangling bond state as well as quasilocalized defect state [4,13]. Previous theoretical studies 
of magnetism in two-dimensional graphitic systems have focused on the graphene monolayer, 
neglecting the interlayer coupling. In this paper, we investigate, using first-principles calcu- 
lations, graphene bilayers with a monovacancy to understand the influence of the adjacent 
graphene layer to the spin magnetic moment of the graphene layer with a monovacancy. 
This result would provide a more realistic understanding of vacancy-induced magnetism in 
graphene multilayer systems including highly oriented pyrolytic graphite (HOPG) [14,15]. 



II. COMPUTATIONAL DETAILS 



We perform ab initio calculations based on the density functional theory with spin po- 
larization. The wave functions are expanded in the double-^ basis set implemented in the 
SIESTA code [16]. Norm- conserving TrouUier-Martins pseudopotentials are employed [17]. 
For the exchange-correlation term, we employ the Cerperley-Aldcr type [18] local density 
approximation with spin polarization (LSDA). 3x3x1 Monkhorst-Pack grids with re- 
spect to 1 X 1 graphene cell are used to sample the Brillouin zone and an energy cutoff for 
real space mesh points is 200 Ry. All coordinates are fully relaxed until the forces of each 
atom are smaller than 40 meV/A. For the Bernal-stacked (AB-type) [19] graphene bilayer 
(composed of 255 carbon atoms) with a monovacancy, two configurations labeled Bi and 
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B2 are considered in the supercell of 19.68 x 17.04 (A^) with 22 A thickness of vacuum as 
shown in figure 1. For Bi, a monovacancy in the upper layer faces a C atom in the lower 
layer along z direction. For B2, on the other hand, a monovacancy in the upper layer does 
not face any C atom in the lower layer. The C-C bond length in one graphene sheet is 1.42 
A, and the interlayer distance is ~3.35 A. With the intention of comparison, a graphene 
monolayer with a monovacancy and 127 C atoms (labeled M) is also calculated. To calculate 
spin and charge distributions of all carbon atoms in the systmes, we carry out the MuUiken 
population analysis [20]. 



III. RESULTS AND DISCUSSION 



Each upper layer in the relaxed bilayer systems (Bi and B2) has nearly the same atomic 
arrangement as M in figure 1. The sheets are almost planar and the atoms near a mono- 
vacancy are displaced owing to the Jahn-Teller effect; for Bi and B2, two of the three C 
atoms around the vacancy are rebonded with bond lengths of 1.78 and 1.79 A, respectively. 
The other C atom with a dangling bond protrudes slightly toward the lower layer. The 
unsaturated danghng bond in B2 is somewhat closer to C atoms in the lower layer than 
that in Bi, so that the perturbation from the interlayer interaction may result in the total 
energy of B2 lower than that of Bi by 0.03 eV, owing to the weakening of the dangling bond 
[5]. Spin (magnetic moment) densities in both of Bi and B2 configurations show notice- 
able differences from that in M. Figure 2 shows perspective view images of isovalue density 
surface plots of the spin density distributions (p-|-(r) — P|(r)) in Bi, B2, and M. Red and 
blue surfaces correspond to spin densities of -|- and — 0.01 e /A^, respectively. For all three 
configurations, figure 2 shows that spin densities are spatially localized at the vicinity of 
the vacancy site and their distributions have the mirror symmetry with respect to y-axis. 
However, spin density satellites to the vacancy sites in Bi and B2 become smaller than those 
in M. Magnitudes of the magnetic moment in Bi, B2, and M manifest this difference more 
quantitatively. Magnetic moments in the upper and lower layers of Bi, obtained from the 
MuUiken population analysis [20], are 1.31, 0.04 ^5, and those of B2 are 1.35 and 0.04 /x^, 
respectively. In comparison, the magnetic moment in M is 1.52 /i^. It means that the 
magnetic moments in the upper layer of Bi and B2 decrease by 14 and 11 %, respectively, 
compared to those in M. 
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To understand the origin of the reduction of the magnetic moment, the interlayer charge 
transfer was checked. We calculated increases in the up-spin, the down-spin, and the total 
charge in the upper layers and their decreases in the lower layers of bilayer systems as listed 
in table 1; the charge increase in the upper layers is calculated by subtracting the total 
charge in M from that in the upper layers, and the charge decrease in the lower layers is 
calculated by subtracting the charge in the lower layers from the that of the ideal graphene 
monolayer with 128 C atoms. For the two configurations Bi and B2, 0.18 e and 0.16 e of 
the total charges are transferred from the lower layer to the upper layer, respectively, and 
most of them occupy only energy levels of down-spin density in the upper layer as shown 
in tabic 1. These values of the interlayer charge transfer are associated with the decrease 
of the magnetic moment in the upper layer of Bi and B2, 0.21 and 0.17 /jb, so that we can 
conclude that the reduction of the magnetic moment originates mainly from the interlayer 
charge transfer from the lower layer to the upper layer and their occupation of energy levels 
of down-spin electrons in the upper layer. For the perturbation, the missing C atom at the 
vacancy site is much more important than the C atom with a dangling bond. 2pz orbitals in 
upper and lower layers overlap and are associated with the interlayer interaction (coupling). 
Because of the missing C atom at the vacancy, the local change in the interlayer interaction 
occurs. For Bi configuration, the vacancy site faces a C atom in the other graphene layer. 
On the other hand, for B2 configuration, the vacancy site faces the center of a hexagon in the 
other graphene layer. Therefore, the charge transfer and the magnetic moment difference in 
Bi are somewhat larger than those in B2. 

Next, we study the influence of the interlayer charge transfer on two origins of 
monovacancy-induced magnetic moments: (1) the localized sp^ dangling bond state from 
broken a bonds and (2) the quasilocalized defect state from broken n bonds. To resolve 
them, we analyze contributions of the orbitals (s, p-^, p^^, and p^) in the C atoms in Bi, B2 
and M to the magnetic moments. In figure 3(a), it is shown clearly that the difference be- 
tween the magnetic moments of Bi, B2 and M comes from the contribution of the p^ orbital. 
By comparing bilayer systems (Bi and B2) to M, we recognize that the reduction of mag- 
netic moments mostly originates from the quasilocalized defect state. Figure 3(b) depicts 
top view images of the magnetic moment density difference, cross-sectioned at the graphene 
layers with a monovacancy. Black-colored grids represent the upper graphene layers. The 
resultant distribution of spin density differences in either case is not localized at the vacancy 
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site but has a regular triangular pattern, which corresponds to the typical pattern of mag- 
netic moments induced by the quasilocalized defect state [4,21]. Therefore, the interlayer 
charge transfer compensates for the magnetic moment induced by the quasilocalized defect 

state. 

The situation could be more clarified with the density of states (DOS) and the band 
structure as shown in figures. 4 (a)-(c). Figures 4 (a) and (b) show the comparison of band 
structures between Bi and M and between B2 and M, respectively. Left and right columns 
represent the bands of up-spin electrons and down-spin electrons, respectively. The DOS of 
the graphene bilayer shown in figure 4(c) is calculated from the contribution of its upper layer 
only. Red and black colors represent the contribution of up-spin and down-spin electrons 
in the bilayers Bi and B2, and blue and green represent those in M, respectively. For these 
three systems, we find several similar features. The band structures show that localized sp^ 
danghng bond states are ~0.5 eV below the Fermi energy (Ep) and quasilocahzed defect 
states are pinned at Ep [4]. The DOSs reveal that the magnetic moments stem mainly from 
the locahzed sp^ dangling bond state, and the contribution of quasilocahzed defect state 
are minor in all three systems. However, some differences between bilayer systems and M 
are also shown in these DOSs and the band structures. In contrast to M, the interlayer 
coupling enhances the localization character of quasilocalized defect states near Y point 
below Ep, so that more down-spin states are created below Ep (marked by purple dotted 
ellipse), resulting in occupation of transferred electrons from the adjacent layer mostly in the 
down-spin quasilocalized defect state. The DOS of down-spin electrons shows this tendency 
more clearly. The main peak position in the DOS of the down-spin quasilocalized defect 
state in the bilayer systems moved from above Ep to below Ep with respect to M, and the 
quasilocalized defect state below Ep of the bilayer systems has less dispersive band than 
that of M, demonstrating an enhancement of the localization character of the bilayer system. 
The band dispersion in the band structure can show localization characters. It means that 
a flat band corresponds to a localized state. Therefore less dispersive (flatter) bands are 
associated with quasi-localized states with heavy effective masses. Similar to the case of 
down-spin electrons, the band structure of the up-spin quasilocalized defect state in the 
bilayer becomes more localized near Y point, compared to M. However, the band of the 
up-spin quasilocalized defect state near Y is far below Ep, so that there is little increase 
in the DOS of up-spin electrons of quasilocahzed defect state from the interlayer charge 
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transfer. Consequently, enhancement of the locahzation character of quasilocahzed defect 
states near Y point below Ep and occupation of transferred electron from the adjacent layer 
in the down-spin quasilocahzed defect state are the major factors that contribute to the 
reduction of magnetic moments. Other bands in the bilayer system are also affected by the 
interlayer coupling in the graphene bilayer, but their contribution to the magnetic moment 
is negligible. 

IV. CONCLUSION 

In conclusion, we demonstrated, by ab initio calculations, that the interlayer charge 
transfer to down-spin quasilocahzed defect state with the enhanced localization character 
results in the reduction of the magnetic moment of the total magnetic moment in graphene 
bilayers. Our study sheds light on the physical behaviors of sp^-bondcd carbon structures 
with vacancies and may lead to a new avenue of carbon-based spintronics. 
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TABLE I: Up-spin, down-spin and total charge increases in upper layers and their decreases in 
lower layers of Bi and B2. 



Systems 


Char 
total 


ge incres 
up-spin 


ise in the upper layer (e) 
down-spin 


Char 
total 


ge decre 
up-spin 


ase in the lower layer (c) 
down-spin 


Bi 


0.18 


-0.01 


0.19 


0.18 


0.07 


0.11 


B2 


0.16 


0.00 


0.16 


0.16 


0.06 


0.10 
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[Figure Captions 

Figure 1. Relaxed atomic configurations of two kinds of Bernal-stacked graphene bilayers 
with a monovancy labeled Bi and B2, and of graphene monolayers with a monovacancy 
labeled M. Black and gray atoms represent C atoms in upper and lower graphene layers, 
respectively. Insets show the site at which a C atom would be removed to generate a 
monovacancy before the relaxation. 

Figure 2. Perspective view images of isovalue density surface plots of the spin [magnetic 
moment] density distributions for Bi, B2, and M. Red and blue colors correspond to the 
spin [magnetic moment] densities of + and — 0.01 e/A^, respectively. 

Figure 3. (a) Contributions of the orbitals (s, p^;, Py, and p^) in the C atoms in Bi, B2 
and M to the magnetic moments (b) Top view images of the difference in spin densities in 
the graphene layers with a monovacancy. Black colored grids represent atom positions in 
the upper layers of Bi and B2. 

Figure 4. Band structures of Bi, B2, and M in the energy range from Ep—leV to Ep+leV 
in (a) and (b). Inset shows the Brillouin zone of the two-dimensional rectangular-shaped 
superlattice. (c) Density of state of upper layers of Bi and B2, and M. Purple-dotted regions 
represent the contribution of the interlayer charge transfer to the reduction of magnetic 
moment in the upper layers of Bi and B2 with respect to M. 
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Figure 1 
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Figure 2 
S. Choi et al. 
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